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We theoretically investigate strong coupling between a single molecule and a single metallic 
nanoparticle. A theory suited for the quantum-mechanical description of surface plasmon polaritons 
(SPPs) is developed. The coupling between these SPPs and a single molecule, and the modified 
molecular dynamics in presence of the nanoparticle are described within a combined Drude and 
boundary-element-method approach. Our results show that strong coupling is possible for single 
molecules and metallic nanoparticles, and can be observed in fluorescence spectroscopy through the 
splitting of emission peaks. 
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I. INTRODUCTION 

Quantum optics has recently made its way to the 
field of plasmonics.^'^ This is due to the the rapid 
progress in nanofabrication and measurement techniques. 
Recent experiments have demonstrated the controlled 
coupling of single molecules with metallic nanoparti- 
cles (MNPs)^'^'^'^ and metallic surfaces,^ of coupled 
nanoparticles, '''^ and of donor and acceptor molecules ac- 
cross metal films. ^ Possible applications of such hybrid 
molecule- MNP systems range from biosensing-'^^'-'^^ to ac- 
tive plasmonic devices. 

A key element of the quantum- optics toolbox is the 
strong coupling between a quantum emitter and a res- 
onator, where excitation energy is coherently transferred 
between emitter and resonator. Strong coupling was 
first observed for single atoms in high-finesse optical res- 
onators/^'^^ and more recently for various solid state 
systems, such as semiconductor quantum dots^^'-^^ or 
superconductor circuits.^'' Although strong coupling be- 
tween ensembles of molecules, e.g, J- aggregates of dyes, 
with plasmons has been reported, ""^^ it is unclear whether 
the strong coupling regime can be reached for single 
molecules coupled to MNPs. The reason for this lies 
in the intricate interplay of the molecule-MNP coupling 
strength with the molecular relaxation dynamics, which 
becomes heavily altered in the vicinity of the nanoparti- 
cle. 

It is the purpose of this paper to theoretically investi- 
gate the strong coupling regime between a single quan- 
tum emitter, such as a molecule or collodial quantum 
dot, and a single MNP. We start by developing a theory 
suited for the quantum-mechanical description of surface 
plasmon polaritons (SPPs), the coupling between these 
SPPs and single molecules, and the modified molecular 
dynamics in presence of the MNP. We employ a Drude 
framework for the description of the metal dynamics, 
and compute the quantized SPP modes within a bound- 
ary element method approach. ^^'^^ Our results show that 
strong coupling is possible for molecules and MNPs and 
could be observed in fluorescence spectroscopy through 
the splitting of emission peaks. 

We have organized this paper as follows. In section II 



we show how to compute surface plasmon modes within 
a boundary-element-method approach, and introduce a 
suitable quantization scheme for the surface plasmons. 
We also present details of the theoretical description of 
the coupled molecule-MNP system in presence of scat- 
terings. Sec. Ill presents results of our model calcula- 
tions. We explore the strong-coupling regime for a single 
molecule coupled to a MNP, and identify the pertinent 
parameters for strong coupling. We also discuss limita- 
tions of our model. Finally, in Sec. IV we summarize and 
draw some conclusions. 



II. THEORY 
A. Plasmon quantization 

Although SPPs are generally considered as bosonic 
quasiparticles, most theoretical work does not explicitly 
rely on such description. In linear response one can em- 
ploy the fluctuation-dissipation theorem to relate the di- 
electric response to the dyadic Green tensor of Maxwell's 
theory,^^'^^ where all details of the metal dynamics are 
embodied in the dielectric function, which can be ob- 
tained from either experiments^ or first principles calcu- 
lations. This approach is no longer applicable in non- 
linear response. Also the neglect of plasmon relaxation 
at small timescales, as proposed in other work^^, is not 
suited for the investigation of strong coupling, which crit- 
ically depends on the relative importance of coupling and 
SPP dephasing. In this work we thus follow the seminal 
work of Ritchie, s^ where the electron dynamics in the 
metal is described within the hydrodynamic model. For 
the transition metals Ag and Au, electrons with particle 
density no are assumed to move freely in a medium with 
background dielectric constant eo, which accounts for the 
screening of (i-band electrons. ^^'^^ 

The energy of a classical electron plasma is the sum of 
kinetic and electrostatic energy^^'^^ 

^=1 j {^o(V^)' + P^}^iV. (1) 

Here p{r) is the charge density displacement from equi- 
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FIG. 1: (Color online) (a) Spectrum of cigar-shaped Ag MNP. 
The height of the particle is 40 nm and the height-to-diameter 
ratio is 5 : 1. The solid and dashed lines show the calculated 
spectra for the dielectric function of Ref. 23 and the Drude 
form,^^ respectively. The triangles at the bottom indicate the 
energies of the plasmon modes (d) and (e), and the arrow at 
~ 1.7 eV indicates the energy of the molecular state 1. (b) 
Non-radiative (solid line) and radiative (dashed-dotted line) 
decay rate for a molecule located at a certain distance from 
the MNP [see panel (c)] in units of the radiative free-space 
decay rate 7°. The inset reports the level scheme used in 
our caluclations. (c) Discretized particle surface as used in 
our calculations, (d-f) Surface charge distribution of SPP 
eigenmodes of lowest energy. Only mode (d) has a nonzero 
dipole moment. 

librium, ^{r) is the electrostatic potential induced by 
/)(r), and ^(r) is the velocity potential, whose derivative 
gives the velocity density v = — V^.^^ Throughout we 
use Gauss and atomic units {e = m = h = 1). For the 
SPPs of our present concern we consider surface charge 
distributions a which are nonzero only at the surface of 
the MNP. As detailed in Appendix A, the Hamilton func- 
tion (1) can be rewritten in a boundary element method 
(BEM) approachi9'29 as 

2no 

Wp [2^(eo + 65) + (eo - 65)^]"' G a} (2) 

Here a is the vector of the surface charges within the dis- 
cretized surface elements (see inset of Fig. 1), G is the 
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FIG. 2: (Color online) Simulation for molecule which is ini- 
tially in state 1 [inset of Fig. 1(b)]. The solid and dashed 
lines show the population of the molecular state 1 and the 
SPP dipole mode [Fig. 1(d)], respectively, for two different 
molecule-MNP distances. We use a molecule dipole moment 
of 10 atomic units, which corresponds to a free-space decay 
time 1/7° of approximately one nanosecond. In weak cou- 
pling (corresponding to the 8 nm distance) the non-radiative 
decay rate can be also estimated from time-dependent pertur- 
bation theory, jnr ~ |^a|^/7o, where gx is the coupling con- 
stant between the molecule and the resonant plasmon mode. 
The dotted line in the figure shows the resulting decay. In 
our simulations we include the 40 plasmon modes of lowest 
energy. 



free Green function matrix which connects two surface 
elements, F is the corresponding surface derivative, ^^'^^ 
ujp = (47rno)2 is the plasma frequency, eo is the back- 
ground dielectric constant of the metal^^, and 65 the di- 
electric constant of the embedding medium. We can now 
determine the eigenmodes of Eq. (2) and quantize the 
plasma oscillations via a canonical transformation, fol- 
lowing the standard procedure outlined in Refs. 22,25,26. 
Within such an approach we obtain the plasmon Hamil- 
tonian H^i = J^a ^a^a second-quantized form, with 
ujx being the energy and a\ the creation operator for the 
plasmon mode A. The field operator for the SPPs is of 
the form 

-W = E(^)'Mr')(a. + 4), (3) 

where ux{r) is the plasmon eigenfunction and /3x the cor- 
responding normalization constant. 



B. Molecule-MNP coupling 

With the SPP quantization we have now opened the 
quantum optics toolbox. This allows us to study strong 
coupling according to the standard prescription^^. As 
for the description of the molecule we follow Refs. 31,32 
who considered a generic few-level system. This approach 
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FIG. 3: (Color online) Fluorescence spectra of the molecule in vicinity of the MNP for different molecular dipole moments, 
given in the panels in atomic units. The moments correspond to free-space decay rates I/7J? of approximately (a) 0.1 /is, (b) 
4 ns, and (c) 1 ns. In the weak-coupling regime of panel (a) the line broadens with decreasing distance but does not split. In 
the strong-coupling regime of panels (b,c) the line splits at a given distance into two lines. 



is also best suited for other quantum emitters, such as 
collodial quantum dots. The inset of Fig. 1(b) shows the 
level scheme used in our calculations. It consists of the 
molecule ground state and two excited states 1 and 
2. We assume that an external pump laser brings the 
molecule into the excited state 2, where it decays non- 
radiativly with a given rate 7^ to the optically active 
state 1. This indirect process allows us to separate the 
excitation dynamics, which is not modified in presence of 
the MNP, from the relaxation dynamics of state 1, which 
becomes strongly modified if the molecule and SPP are 
in resonance. The coherent part of the molecule-MNP 
dynamics is described by the Hamiltonian 



H = H, 



mol 



mol— pi 



pump • 



(4) 



Here i^moi and H^i describe the molecular states and 
the SPP modes, respectively, i^moi-pi is the coupling be- 
tween the molecular dipole and the surface charge (3), 
and i^pump is the interaction of the molecule with the 
pump laser. The last two terms in Eq. (4) are described 
within the usual rotating- wave approximation.^^ In addi- 
tion, we account for the incoherent part of the dynamics 
through a master equation of Lindblad form^^'^^ 



p = -i[H, p]-lYl {l\L^P + ^LiP - 2L^pLl) , (5) 



where p is the density matrix of the coupled molecule- 
SPP system. The Lindblad operators Li describe the 
various scattering channels of molecular decay, plasmon 
decay through Landau damping, and radiative decay^^. 



III. RESULTS 

In our calculations we consider the cigar-shaped Ag 
MNP shown in Fig. 1. Other MNP shapes and metals will 
be discussed at the end. Figure 1(a) shows the spectra 
computed within our BEM approach^^ for the Ag dielec- 
tric function of Ref. 23 (solid line) and the Drude frame- 
work^^ (dashed line). Both spectra are in nice agreement, 
thus justifying the use of the Drude model. The ener- 
gies of the SPP eigenmodes are indicated by triangles. 
Panel (b) reports the non-radiative and radiative decay 
rates of the molecule as a function of molecule-MNP dis- 
tance, which we compute in the weak-coupling regime 
according to the prescription of Refs. 3,20. One observes 
that the rates dramatically increase when the molecule 
approaches the MNP. Here the decay process becomes 
strongly altered by the nanoparticle, which acts as a sup- 
plemental antenna and converts part of the molecule's 
near field into radiation and Ohmic dissipation. 

We next turn to the results of our master-equation 
approach. Figure 2 shows simulations based on the solu- 
tion of Eq. (5) where the molecule is initially brought 
into the excited state 1. Let us first consider the 
larger molecule-MNP distance of 8 nm (upper two lines). 
Through the coupling i^moi-pb the lowest SPP mode be- 
comes populated and subsequently decays through Lan- 
dau damping and radiation. After a transient at early 
times, both molecule and SPP population decay mono- 
exponentially with the same decay constant. In this 
regime the molecule drives the strongly damped plas- 
mon mode and hereby constantly transfers energy to the 
MNP. Things change considerably when the molecule is 
brought closer to the MNP. For z = 4 nm (lower two 
lines) one observes a pronounced population beating be- 
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tween the molecule and the surface plasmon, superim- 
posed on a sub-picosecond decay due to efficient plas- 
mon damping. This beating behavior is a clear signa- 
ture of strong coupling^^ which occurs in a regime where 
the molecule-MNP coupling is stronger than the plasmon 
damping. 

Although strong coupling is most apparent in the time 
domain, spectroscopy appears to be a more suitable tool 
for its experimental observation. We next turn to the 
study of the setup shown in the inset of Fig. 1(b), where 
a weak pump laser brings the molecule to the excited 
state 2. This process is followed by an internal decay into 
the optically active state 1 and a final relaxation to the 
groundstate. Again, the last process is strongly modified 
in presence of the MNP. In our calculations we use the 
master equation (5) to compute the steady state solution. 
Once a stationary condition is reached, we can compute 
the fiuorescence spectra from the Wiener-Khinchin the- 
orem by means of the quantum regression theorem. ^^'^^ 
Figure 3 shows results of our simulations for three differ- 
ent molecular dipole moments. For the smallest dipole 
moment of panel (a), one observes that the line broad- 
ens when the molecule is brought closer to the MNP. We 
verified that the line broadening is precisely given by the 
sum of radiative 7^ and non-radiative 7nr decay rates 
shown in Fig. 1(b). For the larger dipole moments inves- 
tigated in panels (b,c), we observe that at a distance of a 
few nanometers the line splits, thus indicating the onset 
of strong coupling. Here excitation energy is coherently 
transfered between the molecule and the SPP. 

In a generic model, where a quantum emitter is coupled 
to a single cavity mode, the polariton eigenmodes of 
the coupled system are of the form^^ 



^^±+ = ^o-^(7c + 7m)±y5^- (^^^^^)'. (6) 

Here, coq is the enery of the isolated molecule and cavity, 
which are assumed to be in resonance, 7^ and 7c are the 
decay rates of the molecule and cavity, respectively, and g 
is the coupling constant. Strong coupling occurs for g > 
1 7c — 7m 1/4 and corresponds to the formation of a dressed 
state with finite lifetime. It is an intrinsic property of 
the coupling between the molecule and the cavity, and 
manifests itself as a doublet splitting of the emission lines. 
Quite generally, for the coupled molecule-MNP system 
Eq. (6) is too simple, because the molecule couples not 
only to the MNP dipole mode but also to all other modes, 
and one must use a more refined description as we have 
done in this work. Nevertheless, Eq. (6) allows us to 
estimate the pertinent parameters for strong coupling. 
From the plasmon decay rate 70 ~ 30 fs for silver we 
can estimate a critical coupling strength of g ^ jo/A ^ 
5 meV for the onset of strong coupling. Indeed, this 
value is in agreement with the results of Fig. 3 [as can 
be inferred, e.g., from the line broadening in panel (b) 
at the distance of 3 nm where the emission line starts 
to split]. For a coupling of the order of a few meV the 
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FIG. 4: (Color online) Line positions and oscillator strengths 
of polariton modes for cigar-shaped (two lines on left-hand 
side) and disk-shaped (two lines on right-hand side) nanopar- 
ticles and Ag and Au. In Au the lines split at smaller distances 
because of the stronger plasmon damping. We use d = 10 
atomic units. 

approximation of a two- level system is justified for both 
molecules and quantum dots, although the true lineshape 
might be additionally infiuenced by internal degrees of 
freedom (e.g., vibrations) of the quantum emitter. 

When the molecule is brought even closer to the MNP, 
the oscillator strength of the high-energy line vanishes 
and the low-energy line becomes strongly red-shifted. In 
this regime, where the energy renormalization is of the 
order of several tens of meV, the description of the quan- 
tum emitter in terms of a generic few-level system is ex- 
pected to break down. The strong redshift of the emis- 
sion peaks is due to the attractive interaction between 
the molecule and the MNP, which is strongly enhanced 
at small distances, and the different oscillator strengths 
are associated to the different dipole moments of the pre- 
dominantly MNP- and molecule-like polariton modes at 
lower and higher energy, respectively. We also found that 
moderate detunings between the molecule and SPP ener- 
gies do not drastically change the behavior shown in the 
figures. Figure 4 shows that similar behavior is found 
for other MNP shapes and materials. We have also per- 
formed calculations for spherical nanoparticles. Unfor- 
tunately, for both silver and gold the resulting surface 
plasmons have energies in a spectral region where (i-band 
scatterings set in,^'^ and where the Drude description be- 
comes questionable. Our results (not shown) indicate 
that for nanospheres strong coupling occurs at smaller 
distances than for the particle shapes shown in Fig. 4, 
which might be due to the larger number of plasmon 
modes to which the molecule can couple. 

IV. SUMMARY AND CONCLUSIONS 

In conclusion, we have studied strong coupling be- 
tween a single molecule and a metallic nanoparticle 
within a fully quantum-mechanical approach. We have 
demonstrated that strong coupling is possible for real- 
istic molecule and nanoparticle parameters, despite the 
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strong plasmon damping, and should be observable in flu- 
orescence spectroscopy through the splitting of emission 
peaks. Strong coupling is an important ingredient for fu- 
ture plasmonic-based quantum information schemes, and 
might play a significant role in biosensor applications. 
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APPENDIX A 

In this appenix we derive Eq. (2) and show how to 
quantize the SPP modes. Our starting point is the en- 
ergy of a classical plasma, Eq. (1). Let us first consider 
the first term on the right-hand side which describes the 
kinetic energy. From the relation v = — between 
the velocity field v and ^, we obtain for the continuity 
equation 



(Al) 



which gives us the relation between the density displace- 
ment n and the velocity potential ^. In the following we 
consider surface charge distributions a which are nonzero 
only at the surface of the MNP. Integration of the conti- 
nuity equation (Al) over a small cylinder Q (height h ^ 
and base dS)^ which encloses a small surface element, 
then gives for the right-hand side of Eq. (Al) 

/ noV^^(iy= / noii-W'^dS^no^SS. (A2) 
Jn Jon 9n 

Here, ^ = ' denotes the surface derivative of the 
velocity potential. Together with the left-hand side of 
Eq. (Al) we find the link between ^ and cr. 



riQ—r =(j . 
on 



(A3) 



Using Green's first identity we can rewrite the term for 
the kinetic energy in Eq. (1) as 

/ {VmfdV= I ^^dS- I W^^dy. (A4) 

As evident from the continuity equation (A2), for a pure 
surface charge distribution V^^ is zero inside the metallic 
nanoparticle, and the second term on the right-hand side 
of (A4) thus vanishes. We can now use the boundary- 
element method^^'^^'^^ to relate ^ to Our starting 
point is 



^(r) 



/ fG(r, 
Jon \ 



d^s') dG{r,s') \ dS^ 

dh dh ) Att 



where G{r^r') = l/|r — r'\ is the free-space Green func- 
tion. Performing the limit r ^ s in Eq. (A5) according 
to the prescription given in Refs. 19,22,29 and using the 
same notation as in these references, we obtain 



(A6) 



Here and F are the surface derivatives of the velocity 
potential ^ and the Green function G, respectively, and 
G, F and ^, ^' are assumed to be convoluted in space. 

At this point it is convenient to switch to the bound- 
ary elements of the discretized MNP surface (see also 
inset of Fig. 1): ^ and a are vectors of the length of the 
number of surface elements, and G and F are matrices 
connecting the different surface elements. We can thus 
solve Eq. (A6) through inversion ^ = (27r + F)~^G^f' . 
Together with the relation (A3) , the term for the kinetic 
energy can then be brought into the final form 



/ * 



dh 



dS 



^no(7^(27r-hF)-^G(7. (A7) 



Here denotes the transposed surface charge vector. 

For the potential energy of Eq. (1) we follow the pro- 
cedure given in Ref. 19. We start from a relation similar 
to Eq. (A5) but the velocity potential ^ replaced by the 
electrostatic potential ^. Taking the surface derivative 
inside and outside the MNP, we obtain 

(27r-hF)^ = G^[ 
(27r-F)^ = -G$'2- 

Here ^[ and $2 denote the surface derivatives of the 
potential inside and outside the MNP. Multiplying the 
first equation with the dielectric constant eo of the metal 
and the second one with the dielectric constant of the 
embedding medium, gives after substraction 

(27r(eo + eb) -h (eg - eb)F) ^ = G{eo^[ - e^^) = 47r Ga . 

(A8) 

To arrive at the last term of (A8) we have used the 
boundary condition n- {D2— Di) = — e5$2 + ^o^i = 47rcr 
of Maxwell's equation. Putting all above results together, 
we finally get Eq. (2). 

To obtain the eigenmodes we first rewrite Eq. (2) in 
the short-hand notation 



H=--- (a^Ba^a^Aa) , 
2no ^ ^ 



(A9) 



where the explicit form of the matrices A and B can 
be inferred from Eq. (2). The matrices A and B are 
symmetric and thus can be diagonalized simultaneously. 
Let and ux denote the eigenvalues and eigenvectors 
of the generalized eigenvalue problem 

Aux=u;lBux. (AlO) 

The eigenvectors ux can be chosen real, and are orthog- 
onal in the sense 



(A5) 



ul B ux> 



(All) 
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We can thus expand the surface charge distribution in 
terms of these eigenfunctions viz. 

a = ^^xe'^^'axux. (A12) 

A 

Here 7a = y/2no/{uJxPx), and ax is an expansion co- 
efficient for the eigenmode A. Inserting this expression 
into equation (2) and performing the standard quanti- 
zation procedure via a canonical transformation, ^^'^^'^^ 



then brings us to the plasmon Hamiltonian in second- 
quantized form 

A 

with aj^ being the creation operator for the plasmon mode 
A. 
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